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Abstract: 

Measurements of the optical turbulence profile above Siding Spring Observatory were conducted during 
2005 and 2006. This effort was largely motivated by the need to predict the statistical performance 
of adaptive optics at Siding Spring. The data were collected using a purpose-built instrument based 
on the slope-detection and ranging method (SLODAR) where observations of a bright double star 
are imaged by Shack-Hartmann taken with the Australian National University 24 inch and 40 inch 
telescopes. The analysis of the data yielded a model consisting of a handful of statistically prominent 
thin layers that are statistically separated into the ground layer (37.5, 250 m) and the free atmosphere 
(1, 3, 6, 9, 13.5 km) for good (25%), typical (50%) and bad (25%) observing conditions. We found that 
ground-layer turbulence dominates the turbulence profile with up to 80% of the integrated turbulence 
below 500 m. The turbulence tends to be non-Kolmogorov, especially for the ground-layer with a 
power law index of /3 ~ 10/3. The mirror/dome seeing can be a significant fraction of the ground-layer 
turbulence. The median atmospheric seeing, is around 1.2", in agreement with observational reports. 

Keywords: site testing — atmospheric effects — instrumentation: miscellaneous — instrumentation: 
adaptive optics 



1 Introduction 

Of interest is the performance of Adaptive Optics (AO) 
at Siding Spring Observatory (SSO), Australia. The 
motivation for turbulence profiling at SSO is to un- 
derstand the structure of atmospheric turbulence at 
a moderate quality astronomical site and to determine 
the performance predictions for AO. It is know that the 
performance of AO strongly d epends on th e structure 
of the atmospheric turbulence ( Hardy|19 98). No previ- 
ous detailed site-testing of the structure of atmospheric 
turbulence (strength and speed) has been undertaken 
at SSO. The improvements in cost, availability and 
technology make AO a worthy study at such moder- 
ate seeing sites like SSO. AO significantly improves 
the image quality by compensating for the aberrations 
induced by atmospheric turbulence in real-time. The 
performance of an AO system can be predicted with 
simulation codes using input atmospheric model opti- 
cal turbulence profiles (model-OTP) that characterise 
the turbulence above the astronomical site. It could be 
that the installation of adaptive optics for the 3.9 m 
AAT may open the door for new science programs and 
discoveries that would lead to better science. 

If the bulk of the turbulence is low at SSO, then 
Ground-Layer Adaptive Optics (GLAO) correction mode 
can provide significant gains. GLAO provides a large 
FoV (6 arcminutes), but with a partial AO correc- 
tion ( Hubin et al.|2006 i. Certain science cases (includ- 
ing cosmology and extra-galactic observations) require 



larger FoVs with excellent seeing conditions, which 
can be achieved wi th GLAO for a larger fraction of 
nights (iHubin et all 120061). The GLAO correction 



mode averages the wavefront from several or more widely 
separated wavefront sources and applies the result to 
a single deformable mirror conjugated at the ground. 
Hence GLAO performance is typically best when the 
bulk of turbulence is near the ground. It is known that 
the bulk fraction of atmospheric turbulence at most 
astro nomical sites is located near the ground ( |Hardy 
19981, in the "boundary layer" below 500 m in alti- 
tude. The figure of merit for the GLAO correction 
mode is the Ensquared Energy (EE) fraction in a pixel. 
In some cases the improvement in Ensquared Energy 
can be more than double the natural seeing and hence 
halves the integ ration times to ach ieve the same signal- 
to- noise ratios ( Hubin et al.|2006 l. 

Hence, the observations of turbulence profiles are 
of critical importance as the statistical analysis reveals 
a set of model optical turbulence profiles that serve as 
input atmospheres into AO simulation codes. This fact 
has provided the necessary justification for turbulence 
ranging campaigns at major astronomical sites as well 
as research into various site-testing instruments. 

The characteristic structure of the atmospheric tur- 
bulence above SSO was not well understood prior to 
our site-testing campaign. However, the seeing result- 
ing from the total turbulence integral is better under- 
stood wit h DIMM seei n g mea surements at SSO re- 
ported by Wood et al. ( |1995 I. The common seeing 
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values reported by |Wood et al. (19951 are 1.2". Sev- 
eral turbulence profiles have been observed above SSO 
(January 1997) using the Generalized SCIDAR method 
with the ANU 40 " telescope and are reported by |Kluec^ - 
ers et al. ( 1998 1 where is is noted that the strongest 



turbulent layers are located near the ground in the so 
called boundary-layer with heights below 3 km. How- 
ever, a statistically robust model-OTP cannot be de- 
riv ed due to an insufficien t number of profiles observed 
by |Klueckers et al.|p998| . 

This paper discusses the C%(h) and V(h) profile 
measurements taken at SSO and the derived model- 
OTP for suitable use in AO simulation codes. Section 
2 introduces the turbulence parameters. Section 3 out- 
lines the technique and instrumentation used in the 
measurement of the C%(h) and V(h) profiles. Section 
4 and 5 discusses the data and trends noted in the mea- 
sured profiles. Section 6 introduces the model-OTP 
that characterises the atmospheric turbulence profile 
at SSO. Concluding remarks are in section 7. 

The predicted performance of AO at SSO based 
on the presented model-OTP will be published in a 
forthcoming paper. 



2 Turbulence Parameters 

Atmospheric turbulence exhibits a physical process that 
is complex and random in nature, requiring a suit- 
able model. A wi dely accepted model is that proposed 
by Kolmogorov (Hardy 19981, who investigated the 



mechanical structure of atmospheric turbulence. The 
Kolmogorov model described the velocity of motion in 
a fluid medium, where energy is added in the form 
of large-scale disturbances which then break down to 
smaller and smaller structures, until an inner-scale is 
reached. 

From the Kolmogorov model, a spatial power spec- 
trum of phase (power law index, /3 — 11/3) can be 
deduced followed by a set of structure functions, that 
describe non-stationary random fluctuations encoun- 
tered by atmospheric turbulence. Using these struc- 
ture functions, a set of general turbulence parameters 
can be specified that summarize the effects of atmo- 
spheric turbulence. The key parameter used in the 
calculation of the turbulence parameters is the refrac- 
tive index structure constant, C% (units m~ 2 ^ 3 ), and 
its variation with altitude, z, and time, t. We now 
proceed with the specification of the most useful tur- 
bulence parameters for AO. 

The full width half maximum (FWHM) seeing an- 
gle or image dispersion for long exposure images, 4> 
(rad), is given by: 



ro 



(1) 



where A (m) i s the wavele ngth and ro (m) is the 
coherence length ( Fried|1966 1 . Partial AO compensa- 
tion results in a diffraction limited core surrounded by 
a broader halo equal to the seeing disk, <f). The param- 
eter ro includes the integrated effect of the refractive 
index fluctuations over the vertical propagation path, 
z (m). It represents a fictitious "cell size" of turbu- 



lence and defines an aperture diameter over which the 
mean-square wavefront error is 1 rad 2 given by: 



r = 0.185A /5 (sec<C)" 



C n (z)dz 



(2) 



where ( (rad) is the zenith angle. The coherence 
length, ro, also corresponds to the approximate spatial 
scale that adaptive optics must measure and compen- 
sate the effects of atmospheric turbulence. The an- 
gular displacement over which the mean-square wave- 
front error is 1 rad 2 is called the isoplanatic angle 
(Fried||1982[), 6 (rad), given by: 



9 = 0.058A /5 (secC)~ 



Cl(z)z' k {z)dz 



(3) 

The isoplanatic angle, 9q, can be considered the ap- 
proximate compensated field of view for on-axis single- 
conjugate adaptive optics; or the maximum allowable 
angular distance from the high-order wave front source 
to the science object. The coherence time ( Greenwood 
|1977[ ), ro (s), is roughly the time taken for the wind to 
move turbulence by ro. This is given by: 



r = 0.058A % (secC)~ 



Cl{z)vl\ nd {z)dz 



(4) 

The coherence time, ro, can be considered as the 
maximum duration that the atmospheric turbulence 
can be considered 'frozen'; or the maximum duration 
allowable between sequential wavefront samples and 
corrections for the adaptive optics control system. 

Common site-testing instruments that measure at- 
mospheric turbulence usually assume a Kolmogorov 
model of turbulence with power law index, /3 = 11/3. 
All of these parameters can also be measured in the 
case of non-Kolmogorov turbulence but they all be- 
come functions of /3. To estimate the performance 
of optical systems in non-Kolmogorov t urbulence, the 
power spectral density can be expressed ( Stribling et al. 
19951) as 



=a(P)B(z)K-P, (5) 

where Q n (n, j3, z) is the power spectral density as a 
function of position, z is along the optical path, /3 is 
the power law slope (11/3 for Kolmogorov), B(z) is 
the index structure constant having units m 3- ^ and 
a(/3) is a function to maintain consistency with the 
index structure function and is given by Stribling et al. 
(19951 as 



D n (r) = B( z y- 3 . (6) 

The analysis of the non-Kolmogorov model is ap- 
propriate given that the value of the power law index, 
P, can be determined by the fun ction fitting SLODAR 
method ( |Butterley et al.||2006 l. 
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3 Turbulence Measurement 



Various methods are used for turbulence profiling, in 
eludin g direct sensing with microt hermal sensors on 
towers 



Pant et al. (f999) or balloons 



Azouit and Vernin 



(20051), remote-sensi ng wit h acoustic scattering (SO- 
LI AR) Travouillon (|2006 1 , or triangulation of scin- 
tillation (SC IDAR) qVernin and Roddier||T973| |Fuchs 
et al.|[l994t or of image motion (SLODAR) ||Wil- 
son|2002| |Butterley et al.|2006| [Goodwin et al.|2007 1 



These techniques have reached a degree of maturity ex- 
hibiting rea sonable agreement when us e d toge ther in 
campaigns ( Tokovinin and Travouillon] ( 2006 1 , Cerro 
Tololo campaign (ISarazin et al. 2005 I). Each technique 



has its unique benefits and limitations in terms of cost, 
height resolution, height range, temporal resolution, 
ease of implementation and data reduction complex- 
ity. 

3.1 SLODAR Method 

The SLODAR (Slope Detection And Ranging) tech- 
nique has been used on large telescopes at the ORM, 
La Palma and later on a portable, stand alone, turbu- 
lence profiler for ESO, based on a 40 cm telescope with 
an 8 x 8 Shack-Hartmann Wavefront Sens or (SHWFS) 
(5 cm sized sub-apertures) (|Wilson 2002]). The larger 
apertures 5 to 15 cm of SLODAR relaxes exposures 
times to 4 to 8 ms (typical wind crossing timescales) 
providing more suitable observational targets. The 
ground-layer can be measured with sufficient height 
resolutions (50 to 100 m) by observing widely sepa- 
rated double stars (Wilson 2002) whereas higher alti- 
tudes can be investigated by observing more narrowly 
separated doubles. For these reasons, we have selected 
the SLODAR method for our site-testing campaign 
to measure and characterise the turbulence profiles at 
SSO. 

The SLODAR method, illustrated in Figure [TJ is 
an optical triangulation method with the turbulence 
information extracted from the cross-covariance of the 
wavefront slopes of a double star measured using a 
Shack-Hartmann Wavefront Sensor (SHWFS). The two 
optical wavefront components of the double star with 
separation 6, passing through a single turbulent layer, 
at altitude H, produce copies of the aberrations at the 
telescope pupil that are displaced by S = H8 along the 
axis of the double star separation. The corresponding 
turbulent layer shows up in the spatial cross-covariance 
of the optical wavefronts at a spatial offset 5*. 

The SHWFS measures the wavefront slopes by op- 
tically dividing the telescope pupil into an (TV x TV) 
array of square sub-apertures, accompanying each lens 
in a microlens array and measuring the centroids of the 
spot displacements, being proportional to the averaged 
wavefront slope. The sub-apertures and the detector 
have a sufficient field of view to measure the (TV x TV) 
array of spot patterns from both components of the 
double star simultaneously. The exposure times are 
typically 4- 8 ms to freeze the turbulence, being di- 
rectly proportional to sub-aperture size, w, related to 
the wind speed, v, with crossing timescales, r = w/v. 
The sub-aperture sizes are designed to be approxi- 




H„ =D / 9 



5h=w / 9 



Telescope Pupil 



Figure 1 Diagram illustrating the geometry of the 
SLODAR method for a N=4 system. 8 is the dou- 
ble star angular separation. D is the diameter of 
the telescope pupil and w the width of the sub- 
aperture of the Shack-Hartmann Wavefront Sensor 
(SHWFS) array. The centers of the altitude bins 
are given by ASh where A is the lateral pupil sep- 
aration (units of w) and Sh = w/9. The ground- 
layer can be analyzed in higher-resolution by uti- 
lization of double stars having larger 8. 



mately equal to or less than ro, or ranging from 5 cm 
(poor seeing) to 15 cm (good seeing) depending on the 
median seeing. 

The height resolution is uniform, given by Sh — 
w/9 (at zenith). The highest sampled layer, hn-i = 
(TV — l)8h w H m ax = D/9, where TV is the number 
of sub-apertures across the telescope pupil, with the 
ground layer denoted as ho — with resolution Sh/2. 
The vertical resolution and maximum sample height 
are scaled by the inverse of the air mass, \, or cos(£), 
where £ is the zenith distance. 

3.2 SLODAR Instrument 

The SLODAR site-testing campaign to characterise 
the atmospheric turbulence above SSO consists of re- 
sults that were obtained from 8 one-week observing 
runs spanning years 2005 to 2006 with the purpose- 
built 7x7 (Runs 1-6) and 17x17 (Runs 7-8) SLODAR 
instrument configurations using the ANU 24" and 40" 
telescopes. 

The SLODAR instrument can be compared to that 
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of a SHWFS as used in adaptive optics to measure the 
aberrated optical wavefront. However, the SHWFS of 
the SLODAR instrument has a much wider field of 
view requiring the simultaneous measurement of dou- 
ble star optical wavefront gradients. The double stars 
observed by SLODAR can have angular separations 
up to several arcminutes. A functional diagram of the 
first ANU 7x7 SLODAR instrument installed on the 
ANU 24" telescope at SSO is shown in Figure [2] The 
optical diagram for the 17x17 SLODAR instrument 
on the ANU 40" telescope is similar to that shown in 
Figure [2] except for the 2x optic before the image in- 
tensifier. The 17x17 SLODAR instrument is shown in 
Figure [3] 

An advantage of the SLODAR instrument is a rel- 
atively simple optical design. A single collimator lens 
images the telescope pupil onto the (Microlens Array) 
ML A as shown in Figure [2] The collimating lens per- 
forms the following functions: (i) de-magnifies and col- 
limates the incoming light beams of each double star 
component (i.e. point source at infinity); (ii) re-images 
the telescope pupil onto the MLA. 

The MLA performs the function of optically divid- 
ing the telescope pupil into an array of sub-apertures. 
Each sub-aperture forms two spots (image of double 
star) at their focal plane with image motion caused by 
aberrations in the double star optical wavefronts. The 
width of the sub-apertures, w, are comparable to the 
seeing coherent length, ro, keeping the aberrated opti- 
cal wavefront approximately linear. The sub-aperture 
spot displacement from mean position (centroid) is 
proportional to the averaged wavefront gradient (or 
slope). The images of the telescope pupils formed at 
the MLA are completely overlapped and are produced 
by double star components A and B. The individual 
lenslets (or sub-apertures) typically have slow focal ra- 
tios that image double star A and B components onto 
the image intensifier input (photocathode). 

The image intensifier performs the function of ap- 
plying a high gain (> 1000) to the incoming signal 
photons overcoming the high read noise of the high- 
frame rate cameras used with the SLODAR instru- 
ments. The image intensifier allows fainter double 
stars to be observed, down to limiting magnitudes in 
the V-band of approximately 5 to 6 (compared to 1 
to 2 without the image intensifier). The image inten- 
sifier re-images and de-magnifies the MLA spot pat- 
terns onto the camera detector. Hence the SLODAR 
instrument has two focal planes requiring correct ad- 
justments during the calibration process. 

The camera detectors used with the SLODAR in- 
struments have the key features (i) large format cam- 
eras (e.g. 1018x1008 pixel array) to image wide dou- 
ble stars; (ii) high frame rates (15, 20, 30, 200 fps) 
for layer wind speed measurements; (iii) short cam- 
era exposures of 2 ms to 8 ms in order to 'freeze' the 
turbulence. The camera exposures are on time scales 
equivalent to the turbulent layer wind crossing times 
of the sub-apertures. 

Complementing the 17x 17 instrument was the purpose- 
built real-time software, (a graphical user interface 
(GUI) application) to satisfy the requirements for in- 
strument control, processing and diagnostics and au- 




Figure 3 The ANU 17x17 SLODAR Instrument 
(5.8 cm sub-apertures) on the ANU 40" telescope 
at SSO (Photograph taken 18 June 2006). 

tonomously logging of centroid data. The real-time 
software significantly increased the number of quality 
observed datasets. 

4 Data Description 
4.1 Observational Plan 

The scheduled observing runs for the SSO turbulence 
profiling are listed in Table [T] The observational plan 
was to sample the turbulence profile for each season 
over the course of the whole year. Turbulence profiling 
results were obtained from 8 one-week observational 
runs spanning years 2005 to 2006 with the 7x7 (Runs 
1-6) and 17x17 (Runs 7-8) SLODAR instruments us- 
ing the ANU 24" and 40" telescopes respectively. Sam- 
pling over three consecutive months was conducted 
from November 2005 to January 2006 to examine any 
possible monthly variation in turbulence profiles. The 
first observational run in May 2005 provided a test run 
of the 7x7 SLODAR instrument on the ANU 24" tele- 
scope. A total of six observing runs were conducted 
with the ANU 7x7 SLODAR instrument, the final run 
being in January 2006. The instrument changeover 
for the observational run in April 2006 to the 17x17 
SLODAR instrument on the ANU 40" telescope was a 
response to the scientific requirement of needing more 
height sampling resolution bins and the measurement 
of the turbulent layer wind speeds. 

The observational list of double star targets are 
tabulated in Table [2] The limiting magnitude of the 
instruments, V ~ 5.5, limited the observational list 
to a couple (on occasions only one) of suitable tar- 
gets for any given time at the telescope. The ability 
to alternate between ground-layer and free-atmosphere 
sampling is facillitated by switching between widely to 
narrowly separated double star targets. The increased 
number of profiles (datasets) for the final two observ- 
ing runs was a result of changing from the manual 
process of logging data to automated logging with the 
introduction of the real-time software. The final run 
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Figure 2 The expanded optical block diagram of the 7x7 SLODAR instrument (first version) that 
attaches via bayonet mount to the focus of the ANU 24 inch telescope at SSO. Diagram shows 5x5 
SHWFS for simplification. 



Table 1 S LODAR observational runs at SSO 



Run 


Scheduled Dates 


Nights 


Profiles 


Telescope /Instrument 


1 


2-8 May 2005 


4 of 7 




ANU 24" 


/ 7x7 SLODAR 


2 


14-20 June 2005 


3 of 7 


31 


ANU 24" 


/ 7x7 SLODAR 


3 


21-27 September 2005 


3 of 7 


86 


ANU 24" 


/ 7x7 SLODAR 


4 


1-7 November 2005 


1 of 7 


16 


ANU 24" 


/ 7x7 SLODAR 


5 


12-18 December 2005 


2 of 7 


37 


ANU 24" 


/ 7x7 SLODAR 


6 


18-24 January 2006 


5 of 7 


136 


ANU 24" 


/ 7x7 SLODAR 


7 


11-17 April 2006 


5 of 7 


450 


ANU 40" 


/ 17x17 SLODAR 


8 


15-21 June 2006 


6 of 7 


1892 


ANU 40" 


/ 17x17 SLODAR 



target (or 'group' datasets having similar time-stamps 
and height sampling) during the seventh observing run 
(11-17 April 2006) and eighth observing run (15-21 
June 2006). Note we use the convention that h=0 km, 
defines the height of telescope primary mirror. These 
observing runs had the highest number of datasets 
logged (see Table [l]) and provide a useful visual indica- 
as shown in Figure|4|and are processed either by MATLAB®tor of the spatial-temporal evolution of the turbulence. 



in June 2006 witnessed the full operation of the real- 
time software which captured 1892 datasets for off-line 
processing. 

4.2 Data Acquisition 

The SLODAR instrument delivers raw camera frames 



( |MathWorks"||2005| (ANU 7x7 SLODAR instrument) 
or by the SLODAR real-time software (ANU 17x17 
SLODAR instrument). The SLODAR real-time soft- 
ware has the capability to log raw camera frames but 
typically only the centroid data is logged. 

A summary of the observational data statistics are 
listed in Table [I] Data is acquired at 15 fps (Pulnix® 
( JAI Inc|2009l> TM1020, 1018 x 1008 pixels) and 30 fps 
(Pixelink® ( |PixeLINK||2009| ) A741, 1280x1024 pix- 
els) to image larger separated double stars and 200 fps 
(Pulnix TM6740GE, 640x480 pixels) for high tempo- 
ral sampling. A typical dataset consists of a minimum 
of 600 camera frames (15 fps) or 4000 camera frames 
(200 fps) for sufficient statistical sampling of the at- 
mospheric turbulence. 

5 Data Analysis 
5.1 Cf](h) profiles 

In this section we provide examples of consecutive C%(h) 
turbulence profiles taken from the same double star 



Consecutive profiles for the seventh run are plotted in 
Figure [5] Example individual turbulence profiles from 
these group datasets are plotted in Figure [6] 

The temporal plots of the turbulence profiles shows 
a number of dynamical characteristics; (i) intense tur- 
bulence occurring near the ground (below 100 m); (ii) 
turbulent layers that can fluctuate in intensity, ap- 
pearing in 'bursts' with timescales of several minutes. 
This can have implications for different lines of sight; 
(iii) appearance to drift in altitude on some occasions 
rather than disperse before disappearing. 

To quantify the structural distribution of the atmo- 
spheric turbulence two statistical parameters are de- 
fined, namely h nnnn . The h nnnn parameter describes 
the fractional amount (0.0 to 1.0, with 1.0 being 100%) 
of turbulence below nnnn m as measured from the tele- 
scope primary mirror. The default value for nnnn is 
500 m or h^oo, but due to coarse height resolution sam- 
pling, the value of nnnn can be larger, e.g, 750 m or 
/1750. The h nnnn allows a qualitative assessment for 
the performance of ground layer adaptive optics which 
is favorable if the bulk of the turbulence is near the 
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Table 2 Double star targets for SLODAR observations at SSO. 



Name 


a (RA) 


5 (DEC) 


mag 


separation 


6 Eri 


02 58 


-40 18 


3.4/4.5 


8.2 " 


DUN 16 Eri 


03 49 


-37 37 


4.8/5.3 


7.9 " 


9 Ori 


05 35 


-05 25 


4.9/5 


135 " 


/3 Mori 


06 28 


-07 02 


4.7/5.2 


7.2 " 


oc Cru 


12 26 


-63 06 


1.3/1.8 


4.3 " 


jU Cru 


12 55 


-57 11 


4.0/3.5 


34.9 " 


HIP65271 


13 23 


-60 59 


4.5/6.1 


60 " 


a Cen 


14 39 


-60 50 


-0.01/1.33 


9.5 " 


9 Ser 


18 56 


+04 12 


4.6/4.9 


21.6 " 


<5 Aps 


16 20 


-78 42 


4.7/5.3 


102.9 " 




Figure 5 SSO Run 7: Examples of consecutive SLODAR turbulence profiles from difference spaced double 
star targets (height ranges) . Each temporal plot represents a group of datasets from the same double star 
target (similar height sampling) measured during 11-17 April 2006. The vertical axis denotes height (km) 
and the horizontal axis denotes time (UTC). The color denotes turbulence strength, Cjj(h).dh (m 4 ~P). 
The plots are derived by interpolating the turbulence profiles onto a regular spaced grid at approximately 
Nyquist sampling. The blank regions represent times having no data. Note, h—0 km, defines the height 
of telescope primary mirror. 
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Figure 6 SSO Run 7: Examples of individual SLODAR turbulence profiles with each plot representing a 
single profile as represented in temporal plots of Figure [5] measured during 11-17 April 2006. The vertical 
axis denotes turbulence strength, C%{h).dh (m 4-/3 ), and the horizontal axis denotes height (km). Note, 
h=0 km, defines the height of telescope primary mirror. 
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(b) 



Figure 4 An ensemble average over raw camera 
frames of the (a) SLODAR 7x7 Instrument on 
ANU 24" telescope at SSO (b) SLODAR 17x17 
Instrument on ANU 40" telescope at SSO. The 
double stars observed are (a) a Crux and (b) 
a Cen are aligned along the SHWFS x-direction. 



ground. 

The third statistical parameter characterises the 
power-law slope index of the power spectrum of spa- 
tial phase fluctuations, f3 a v 3 - The parameter f3 avg is 
representative of the entire atmospheric turbulence, or 
averaged contribution of all layer heights. The /3 avg 
was determined by the best fit of the A = theoret- 
ical covariance impulse response function for ft rang- 
ing from 19/ 6 to 23/6 to the observ ed auto-covariance 
function, see |Butterley et al~ (20061. The implications 
of non-Kolmogorov turbulence (/? ^ H/3) in the mea- 
surement of atmospheric turbulence and astronomical 
imaging are discussed by |Stribling et al.| ( |1995| ) and 
Goodwin|p009 i. 



Fits of a non-Kolmogorov exponent to the SLO- 
DAR power spectrum (auto-covariance) are not pro- 
vided in this paper. The fitted data has comparable 
error bars to |Butterley et al.| ( |2006[ ) and that the best 
fit within the error bars typically produces exponent 
values that are non-Kolmogorov (less than 11/3). An 
exponent less than 11/3 causes the DIMM seeing to 
be overestimated due to increased image motion for 
small apertures |Goodwin| ([2009). The power spec- 
trum is relatively insensitive to the o uter-scale (Von 



Karman power spectrum) as noted by putterley et al 
( |2006[ ) . A fit of the outer-scale would typical need an 
outer scale smaller than the telescope (lm) which is 
not typical of other measurements at other observa- 
tories (20-40m). The outer-scale fit also has typically 
larger residuals for the larger offsets in the covariance 
function compared to the exponent fit. 

The results measured during 11-17 April 2006 (run 
7) and 15-21 June 2006 (run 8) using the ANU 17x17 
SLODAR instrument on the ANU 40" telescope indi- 
cate an atmospheric turbulence structure that is dom- 
inated by strong ground-layer turbulence. This is evi- 
dent in the summary /1500 parameter for all nights as 
shown in Figure [7] The summary /1500 parameter in- 
dicates that nearly 76% (run 7) and 91% (run 8) of 
the integrated turbulence is below 500 m. Note that 
a few cases of limited height resolution sampling, the 
ftsoo will be more representative of higher altitudes, 
e.g. /1750, but this is a minority of the datasets. 

The summary /3 aV g parameter for all nights is shown 
in Figure [8] The summary /3 avg parameter is found to 
have a median of 3.32 (run 7) and 3.35 (run 8) as com- 
pared to a Kolmogorov value of 3.67. This implies that 
the strong ground-layer is non-Kolmogorov causing a 

low fiavg- 

The summary seeing (for a wavelength of 0.5 mi- 
crons) derived from the SLODAR function-fitting method 
(non-Kolmogorov analysis) and the DIMM method (Kol- 
mogorov) are shown as a histograms in Figure [9] The 
median seeing values for non-Kolmogorov and DIMM 
analysis (using SLODAR data) are 0.77" and 1.13" 
(run 7) and 1.1" and 1.33" (run 8). The discrep- 
ancy between SLODAR (non-Kolmogorov) and DIMM 
(Kolmogorov) seeing calculation methods are most likely 
due to the low f3 avg values in the data (strong non- 
Kolmogorov effects). It is important to note that the 
seeing va lues have the mirror/ dome seeing component 
removed ( Goodwin et al.|2007 l, which is usually found 
to be a significant component. We note that the DIMM 
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(a) (a) 




(b) (b) 



Figure 7 Summary plot showing the fraction of tur- 
bulence below 500 m based on all observable nights 
measured during (a) SSO Run 7: 11-17 April 2006 
with the median fractional amount of turbulence 
below 500 m is 76% (based on 450 datasets); (b) 
SSO Run 8: 15-21 June 2006 wit the median frac- 
tional amount of turbulence below 500 m is 91% 
(based on 1892 datasets). 



Figure 8 Summary plot showing the average power 
law slope, P aV g, of the spatial power spectrum of 
phase fluctuations based on all observable nights 
measured during (a) SSO Run 7: 11-17 April 2006 
with the median of 3.32 (based on 450 datasets); 
(b) 15-21 June 2006 with the median of 3.35 (based 
on 1892 datasets). For both cases the values are 
noticeably less than the Kolmogorov value of 3.67 
(dashed vertical line). 
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seeing of the SLODAR data, 1.13" and 1.33" brackets 
the historical DIMM seeing measurements by |Wood| 
|et~aTl ( p95| of 1.25". 



5.2 V(h) profiles 

The translational velocity information of the turbulent 
layers can be retrieved by introducing gradual tempo- 
ral offsets, St, in the spatial cross-covariance function 
of wavefront gradients and observing the correspond- 
ing displacement of peaks ( |Wilson | |2002| |. The tem- 
poral offsets, St, are integer multiples of the camera 
acquisition time (inverse of frame rate), and therefore 
must be sufficiently short to capture multiple obser- 
vations of the turbulent layer as it moves across the 
telescope pupil. Typically, St < 100 ms, with cam- 
era acquisition times between 5 ms and 50 ms. A 
turbulence layer moving at velocity, v, will have its 
cross-covariance peak shifted by vSt from its location 
at St = 0, aligned along the separation axis of the dou- 
ble star. By making several measurements of the spa- 
tial shifts in the cross-covariance peak for several se- 
quential temporal offsets, St, it is possible to trace the 
turbulent layer back to the origin to determine both 
height and velocity information. Examples of a layer 
wind speed measurement using the temporal spatial 
cross-covariance of centroid data during the seventh 
observing run (11-17 April 2006) is shown in Figure 10 



SLODAR tredJ=0.7S5. SD1MM £blllQ=1 .13. N= 450 




(a) 



6 Model-OTP 



A model-OTP is required to summarize the main char- 
acteristics of the measured atmospheric turbulence so 
that adaptive optics simulations can be performed and 
instrument performance can be predicted. The diffi- 
culty is that the atmospheric turbulence profiles follow 
a non-stationary process and therefore individual pro- 
files are not representative for use in adaptive optics 
simulations. Of particular interest are the character- 
istics of the ground-layer and free-atmospheric turbu- 
lence, such as the contribution to the total turbulence 
integral (seeing) , thickness and intensity of the ground- 
layer and if any persistent prominent layers exist. 

Model-OTPs to characterise the atmospheric tur- 
bulence have been synthesized from measurements ob- 
tained at other astronomical observatories. Tokovinin 
and Travouillon ( 2006 ) note that previous model-OTPs 



that are based on the averag e or median profiles, such 
Abahamid et al. (2004), do not model the strong 



variab ility property of turbulence. |Abahamid et al] 



(2004) point out that the turbulence intensity at any 



given altitude changes by several orders of magnitude 
and that real OTPs are typically dominated by a few 
strong layers. Hence the use of only median or average 
techniques for OTPs to characterise the atmospheric 
turbulence for adaptive optics analysis may be mis- 
leading. 

An example is the characterisation of the OTP 
above the Cerro Pachon (CP) astronomical site. The 
CP site was initially characterised during the 1998 
Gemini site campaign ( Vernin et al. 1998 1 in which 
seven discrete layers were modeled by Ellcrbroc k and| 



SLODAR (redt-1.1 S DIMM (blue)- 1.33 U- 




(b) 



Figure 9 Summary plot showing the seeing his- 
tograms of the SLODAR (non-Kolmogorov, red) 
method and DIMM (Kolmogorov, blue) method 
based on all observable nights measured during 
(a) SSO Run 7: 11-17 April 2006 (based on 450 
datasets); (b) 15-21 June 2006 (based on 1892 
datasets). For both cases the seeing valuesare re- 
ported at a wavelength of 0.5 microns. 
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St = ms 5t = 5 ms 5t = 1 ms 




Figure 10 SSO Run 7: Example temporal spatial cross-covariance of centroid data for layer wind speed 
determination as measured during 11-17 April 2006 with the CCD camera (TM6740GE) having a frame 
rate of 200 fps (area of interest read-out). The example is of the double star a Cen with height resolution 
5h=1.0 km. The double star separation axis (positive heights) is marked with a black line. The temporal 
offset, r, is a multiple of the inverse of the frame rate, or 5 ms, starting from top left panel with t=0 ms 
and the largest offset located at bottom right panel with t=40 ms. The pixels represent the sub-aperture 
offsets (8i, 5j) with physical size u>=5.8 cm. The wind speed of a layer can be estimated by s/r, where 
s is the physical displacement of covariance peak for a given temporal offset, r. For this example, four 
separate layers are detected with speeds (1) 0.5 m/20 ms = 25 m/s (8 km); (2) 0.4 m/40 ms = 10 m/s 
(4 km); (3) 0.26 m/40 ms = 6.5 m/s (2 km); and (4) 0.18 m/40 ms = 4.5 m/s (0 km). Note, h=0 km, 
defines the height of telescope primary mirror. 
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Rigaut (20001, referred to as the ER2000 model. The compariso ns are possible. The |Tokoviniri and Travouil-| 



ER2000 model has been used by other groups in AO 
simulations but is insufficient in two respects as noted 
by |Tokovinin and Travouillon (20061: (i) it does not 
address the variability of the real OTP and (ii) it 
was developed for the needs of classical and multi- 
conjugate AO that is mostly affected by high-altitude 
turbulence. 

To resolve the limitations with the ER2000 model, 
and o ther similar models, |Tokovinin and T ravouillon 
( 2006 1 propose a new method to derive a more de- 
tailed statistical model-OTP suitable for GLAO anal- 
ysis as well as other adaptive optics techniques. The 
proposed method has already been used by |Ander-| 
sen et al.| ( |2006| ) to model the OTP for Cerro Pachon 
for wide-field GLAO simulation for Gemini-South. An 



Ion 



[2006 1 method calculates the C N turbulence inte- 
gral, J, derived from observations for the ground layer 
(Jgl) and the free- atmosphere (J fa), then categorizes 
into 'good', 'typical and 'bad' conditions. The repre- 
sentative 'good' profile is based on an averaged profile 
representative of the 1 st quartile (25%) of Jgl and 
J fa, based on averaging observational profiles in the 
range (15% to 35%) to ensure an adequate sample size. 
Likewise the 'typical' profile based on the range (40% 
to 60%) and the 'bad' profile (65% to 85%). This pro- 
cess results in a set of three representative profiles for 
each of the GL and FA that have been averaged sepa- 
rately for each group to reveal typical features. 



The Tokovinin and Travouillon (2006) method fits 



exponential equations to model the GL intensity and 
thickness. In this paper, we adopt a slightly different 
approach in that prominent layers are modeled as thin 
discrete layers, as representative of high-resolution tur- 
bulence profiles with micro-thermal balloon measure- 
ments. The thin layers preserve the relative strengths 
Methodology: Layer Strength Models well as the turbulence integrals J G l for the GL and 

Jfa for the FA based on the 'good' (25%), 'typical' 
(50%) and 'bad' (75%) quartiles. 



outline of the methodology for the model-OTP is given 
by Tokovinin and Travouillon (20061. Note that the 
Model-OTP derived in this paper also assumes a Kol- 
mogorov turbulence model. 



6.1 



The mode l-OTP proposed by Tokovinin an d Travouil-| 



Ion ( 2006 1 separates the ground layer (GL) statistics 



from the free-atmosphere (FA) statistics (observed to 
be independent). The GL zone is defined from the tele- 
scope (10-50 m) to some 500 m above the site. The FA 
zone is defined as all turbulence above the GL zone. 

The Tokovinin and Travouillon ( 2006 1 method is 
well suited to the SLODAR method. The SLODAR 
method has the ability to measure the GL and FA 
with sufficient height resolution by choosing appropri- 
ate double star separations. The GL can be measured 
with sufficient height resolution by selecting double 
stars having the widest separation. 

The OTP is universally defined as the dependence 
of the refractive index structure constant C% (mea- 
sured in units m _2//3 ) at altitude h (measured in units 
of meters) above sea-level (in this paper, h=Q km, de- 
fines the height of telescope primary mirror). The tur- 
bulence integral, J, is defined as: 



J ~ Cjv (h).dh 



(7) 



(measured in units m 1 ^ 3 ) is calculated over some 
altitude range. 

The turbulence integral calculated over the entire 
height range covering the atmospheric turbulence (0 to 
20 km), the total turbulence strength, can be expressed 
as the astronomical 'seeing' (units in arcsecs). The 
'seeing' is the spot image full width at half-maximum 
(FWHM) of an unresolved object (e.g., unresolved star). 
The seeing for observations at a wavelength of A = 
500 nm at the zenith (Kolmogorov turbulence model) 
is given by: 

e= [J/(6.8X10- 13 )]° 6 (8) 

As noted by |Tokovinin and Travoui llon (20061 the 
seeing is not additive, hence the preference to use the 
turbulence integrals, J, which are additive and direct 



The Jgl was calculated using the 'seeing' parame- 
ter (converted to J using Equation[8| and hsoo param- 
eter (fractional turbulence below 500 m), using the fact 
that Jgl — hsoo J- The /1500 parameter can be calcu- 
lated for all turbulence profiles as the maximum height 
is greater than 500 m as well as the height resolution 
of the zero height bin is less than 500 m for almost 
all observations. Hence the /1500 parameter is suitable 
for the calculation of Jgl- The Jfa was calculated 
using the fact that Jfa = J — Jgl- For the simplicity 
of calculations, it was assumed that the GL and FA 
are both modeled by a Kolmogorov turbulence model. 
With three turbulence profiles based on thin-layers for 
each of the GL and FA it is possible to construct an 
OTP model having nine possible outcomes with re- 
spective probabilities. 

The models sufficiently cover half of the conditions 
(25% to 75%) expected at the astronomical site within 
which a probability of 25% is assigned to 'good', 50% 
is assigned to 'typical' and 25% is assigned to 'bad' 
for GL and FA turbulence profiles. The total turbu- 
lence profile is a combination of the GL and FA turbu- 
lence profiles with probability equal to the respective 
GL and FA probabilities being multiplied. The assign- 
ment of probabilities to all possible turbulence profile 
outcomes of the model allows a relative importance to 
be assigned. The model does not sufficiently represent 
the extreme cases of very 'good' (0% to 25%) and very 
'bad' (75% to 100%) conditions. It is noted that the 
very 'bad' condition is of most concern as adaptive op- 
tics may not provide sufficient wavefront correction for 
scientific operations. Sufficient representation for the 
very 'good' and very 'bad' conditions can be obtained 
by extrapolation, scaling the layer relative fractional 
amounts by the Jgl and Jfa turbulence integrals, us- 
ing values from the respective cumulative density func- 
tion (CDF) plots. 
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6.1.1 Ground Layer Model-OTP 

To derive the GL model-OTP we consider three model 
turbulence profiles derived from averaging a group of 
observational profiles in intervals centered on the cu- 
mulative density function quartiles for the Jgl turbu- 
lence integral. The 'good' model turbulence profile is 
derived from averaging profiles having maximum sam- 
pling height, Hmax, greater than 500 m but less than 
2000 m (ground-layer sampling) within the 15% to 35% 
interval of the cumulative values of Jgl- Likewise for 
the 'typical' model turbulence profile, within the 40% 
to 60% interval, and 'bad' model turbulence profile, 
within the 65% to 85% interval. The relatively small 
value of Notp for GL analysis reflects the fact that 
only a limited number of datasets suitable for GL sam- 
pling at SSO were observed. The values relating to the 
GL model turbulence profiles are listed in Table [3] 

From Table [3] we see that some intervals have more 
turbulence profiles, due to the Jgl distribution con- 
sisting of a mixture of both GL and FA turbulence 
profiles, with a portion not meeting the maximum sam- 
pling height criterion. The average maximum height 
range, H ma x, is 0.72 km to 1.07 km and the average 
height resolution, 8h, is 72 m to 108 m. 

Figure [TTJ shows the GL model turbulence pro- 
files obtained from averaging GL profiles within cer- 
tain interval ranges of the Jgl distribution to repre- 
sent 'good', 'typical' and 'bad' seeing conditions, as 
summarized in Table [3] 

A thin-layer GL model-OTP is required to be de- 
rived from the continuous GL model-OTP. A thin-layer 
GL model-OTP is an accurate representatio n of high 
resolution OTPs ( Azouit and Vernin 2005 1 , as well 



as compatible for adaptive optics simulations (using 
phase screens to represent thin layers). For the thin- 
layer GL model-OTP we define two layers at heights 
37.5 m and 250 m to best represent the continuous 
GL model-OTP, Figure [TT] To calculate the strengths 
of the layers an integral (lower, upper) centered on the 
layer heights of the continuous GL model-OTP is per- 
formed. To compare how well layer integrals explain 
the model, the total layer integral (37.5 m + 250 m) is 
compared with the total GL integral from m to 500 m 
(calculated from the continuous GL model-OTP). The 
results are listed in Table [4] 

The fractional strengths can be multiplied by Jgl 
according to cumulative level values obtained from its 
corresponding CDF plot. The model is based on the 
cumulative levels 25%, 50% and 75%, but other levels 
can be used for extreme conditions, eg 10%, 50% and 
90%, providing some flexibility for a custom model- 
OTP. 

The turbulence integral strengths of the model thin 
layers are calculated by multiplying the fractional amounts 
with the cumulative levels 25%, 50% and 75% of Jgl- 
The final model thin layers and their turbulence inte- 
gral, with total ground layer seeing, €gl, and model 
probability, are listed in Table [5] 

6.1.2 Free Atmosphere Model-OTP 

To derive the FA model we follow a similar approach 
to that described for the GL model. This involves 





Model Integral, J 


Layer Height (m) 


good typical bad 


37.5 
250 


4.7129 7.3989 9.4363 
0.2485 1.1574 3.8486 


Total (37.5+250) 


4.9614 8.5562 13.2849 




0.8277 1.1478 1.4945 


Probability 


25% 50% 25% 



Table 5 Final turbulence integrals for the GL thin- 



layer model-OTP, J in units of 10 
(Run 1-8: May 2005 to June 2006). 
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deriving three model turbulence profiles from averag- 
ing a group of observational profiles within represen- 
tative intervals of cumulative density function values 
for the Jfa turbulence integral. The 'good' model 
turbulence profile is derived from averaging profiles 
having maximum sampling height, Hmax, greater than 
16000 m but less than 20000 m (free-atmosphere sam- 
pling) within the 15% to 35% interval of the cumulative 
values of Jfa- Likewise for the 'typical' model turbu- 
lence profile, within the 40% to 60% interval, and 'bad' 
model turbulence profile, within the 65% to 85% inter- 
val. The values relating to the FA model turbulence 
profiles are listed in Table [6] 

From Table [6] we see that some intervals have more 
turbulence profiles, due to the Jfa distribution con- 
sisting of a mixture of both GL and FA turbulence pro- 
files, with a portion not meeting the maximum sam- 
pling height criteria. The average maximum height 
range, H ma x, is 17.0 km to 17.3 km and the average 
height resolution, Sh, is 1095 m to 1338 m. 

Figure fl2] shows the FA model turbulence profiles 
obtained from averaging FA profiles within certain in- 
terval ranges of the Jfa distribution to represent 'good', 
'typical' and 'bad' seeing conditions, as summarized in 
Table H 

A thin-layer FA model-OTP is required to be de- 
rived from the continuous FA model-OTP. A thin-layer 
FA model-OTP is an accurate representation of high 



resolution OTPs (Azouit and Vernin|2005 ) as well as 



compatible for adaptive optic simulations (using phase 
screens to represent thin layers). For the thin-layer FA 
model-OTP we define five layers at heights 1000 m, 
3000 m, 6000 m, 9000 m and 13500 m to best represent 
the continuous FA model-OTP, Figure[l2] To calculate 
the strengths of the layers an integral (lower, upper) 
centered on the layer heights of the continuous FA 
model-OTP is performed. To compare how well layer 
integrals explain the model, the total layer integral 
(1000 m+3000 m+6000 m+9000 m+13500 m) is com- 
pared with the total FA integral from 500 m to 16000 m 
(calculated from the continuous FA model-OTP). The 
results are listed in Table [7] 

The FA model-OTP in Table [7] has a gap in bin- 
ning between 10.5km and 12km as well above 15km. 
This is because we rarely saw any evidence for tur- 
bulence at those heights in the continuous FA profile 
of Figure |12| There was also an attempt to keep the 
number of model layers to minimum as well as keeping 



14 



Publications of the Astronomical Society of Australia 



Table 3 Levels of the cumulative distributions of Jgl used in the calculation of a representative ground 
layer profiles, 'good', 'typical' and 'bad' for SSO (Run 1-8: May 2005 to June 2006). 





lower 


upper 


Notp 




Ih 


Jgl 




Good 


15% 


35% 


17 


1.0728 


0.0815 


5.1866 


0.8483 


Typical 


40% 


60% 


10 


0.7211 


0.1083 


8.7999 


1.1662 


Bad 


65% 


85% 


32 


0.7985 


0.0722 


13.5809 


1.5124 




I — , — I — , — I — , — I — I — I — I — I — , — I — , — I — , — I — , — J. — I — I 02 l 1 1 1 1 1 1 1 1 1 1 

5D 1DD 15D 20D 250 300 350 400 450 500 50 100 150 200 250 300 350 400 450 500 

Altitude, m 



(a) (b) 

Figure 11 Continuous GL model-OTP for SSO (Run 1-8: May 2005 to June 2006): squares 'good'; 
diamonds 'typical' and circles 'bad' (a) averaged profiles, error bars are 95% confidence interval and (b) 
corresponding CDF profiles. 





Range (m) 


Integral, J 


Layer Height (m) 


lower upper 


good typical bad 


37.5 
250 


150 
150 350 


4.6104 6.9491 8.3327 
0.2431 1.0870 3.3985 


Total (37.5+250) 




4.8535 8.0361 11.7312 








Integral (0-500) 


500 


5.1866 8.7999 13.5809 


% Explained 




93.5782 91.3203 86.3805 



Table 4 Turbulence integrals for the thin-layer model-OTP for the ground layer, J in units of 10 13 m 1 ^ 3 
for SSO (Run 1-8: May 2005 to June 2006). 



Table 6 Levels of the cumulative distributions of J F a used in the calculation of a representative ground 
layer profiles, 'g ood', 'typical' and 'bad' for SSO (Run 1-8: May 2005 to June 2006). 





lower 


upper 


N tp 




Ih 


Jgl 




Good 


15% 


35% 


147 


17.0076 


1.0951 


0.3702 


0.1719 


Typical 


40% 


60% 


118 


17.1918 


1.3382 


1.1983 


0.3509 


Bad 


65% 


85% 


58 


17.3392 


1.1926 


2.5276 


0.5503 
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Altitude, in 

(a) (b) 

Figure 12 Continuous FA model-OTP for SSO (Run 1-8: May 2005 to June 2006): squares 'good'; 
diamonds 'typical' and circles 'bad' (a) averaged profiles, error bars are 95% confidence interval (b) 
corresponding CDF profiles. 



equal bin widths of 3000m. The gap in binning will 
result in a slight bias to increase the relative strength 
of the FA lower layers causing the isoplanatic angle to 
be slightly larger (providing an optimistic estimate). 
The value of the bias is expected to be small. 

The fractional strengths can be multiplied by Jfa 
according to cumulative levels obtained from its CDF 
plot. The model is based on the cumulative levels 25%, 
50% and 75%, but for other levels extreme conditions 
can be used, e.g. 10%, 50% and 90%, providing some 
flexibility for a custom model-OTP. 

The turbulence integral strengths of the model thin 
layers are calculated by multiplying the fractional amounts 
with the cumulative levels 25%, 50% and 75% of Jfa- 
The final model thin layers and their turbulence inte- 
gral, with the total free atmosphere seeing, cfa, and 
model probability, are listed in Table [8] 

The total integrated turbulence of the thin-layer 
FA model as shown in Table [8] is approximately twice 
that of the total integrated turbulence of binned layer 
continuous FA model as shown in Table [7] A likely 
explanation of the scale factor of ~ 2 may be the re- 
sult that the parameter Jfa has been overestimated 
(turbulence integral in the free-atmosphere computed 
based on total seeing and -ffsoo) compared with the di- 
rect turbulence integrals of the averaged measured pro- 
files. The total strength of the weaker free-atmosphere 
layers (not location) as measured by SLODAR is some- 
what underestimated (low S/N). The factor of ~ 2 
scaling increase makes the free-atmosphere seeing in 
the final model a conservative estimate. 

6.2 Methodology: Layer Wind Speed 
and Direction Model 

To model the turbulent layer speeds a Button wind 
speed profile model is used (see Equation]!)]). To model 
the 'good', 'typical' and 'bad' conditions, three sepa- 
rate Button wind profiles are presented based on data 



found in literature ( Azouit and Vernin 2005 Tokovinin 



et al.|2003||Avila et al.|2003[ ). The wind profiles are not 
based on fits to our data as it is not always possible to 
obtain a sufficient sample of layer wind speeds with our 
SLODAR data due to (i) low camera frame rates; (ii) 
finite camera exposures; (iii) weakness of a layer; and 
(iv) short boiling lifetimes (layer de-correlates rapidly). 
The Button wind profile model is given by 



v(h) = vg + vrexp 



(9) 



where vc denotes the wind velocity at low alti- 
tude, vt denotes the wind velocity at the tropopause, 
Ht denotes the height of tropopause, Lt denotes the 
thickness of the tropopause layer. 

The parameters of the Button wind model for 'good', 
'typical' and 'bad' conditions used in the analysis are 
listed in Table [9] The wind model associates strong 
ground layer wind speeds with stronger free-atmosphere 
wind speeds that have a broader upper atmosphere ex- 



tent (Azouit and Vernin||2005| [Tokovinin et al.||2003 



Avila et al.|2003 l. To simplify the model, the 'good' 
'typical' and 'bad' conditions are represented by the 
same profile and referenced to the ground-layer wind 
direction, set to degrees. The free atmosphere layers 
can travel in a direction perpendicular to the ground- 
layer direction. The wind direction model, tp(h) for 
'good', 'typical' and 'bad' conditions is defined as: 

ip(h) = o(l - exp(~h/b)) (10) 
where a = 101.93 and b — 3255.4 are used. 



6.2.1 Turbulent Layer Wind Model 

To model the turbulent layer wind speed and direc- 
tion for the SSO (Run 1-8: May 2005 to June 2006) 
runs a Button wind speed model together with layer 
height information from the model-OTP is used. The 
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Range (m) 


Integral, J 


Layer Height (m) 


lower upper 


good typical bad 


1000 
3000 
6000 
9000 
13500 


500 1500 
1500 4500 
4500 7500 
7500 10500 
12000 15000 


0.0770 0.2577 0.3675 
0.0100 0.2358 0.6265 
0.0065 0.0499 0.1974 
0.0093 0.0466 0.1222 
0.0209 0.0641 0.1810 


Total (1000+3000+6000+9000+13500) 




0.1237 0.6541 1.4946 








Integral (500-16000) 


500 16000 


0.1409 0.7032 1.6387 


% Explained 




87.7795 93.0185 91.2053 



Table 7 Turbulence integrals for the thin-layer model-OTP for the free atmosphere, J in units of 10 13 to 1 / 3 
for SSO (Run 1-8: May 2005 to June 2006). 





Model Integral, J 


Layer Height (m) 


good 


typical 


bad 


1000 
3000 
6000 
9000 
13500 


0.2413 0.4772 0.6677 
0.0315 0.4368 1.1384 
0.0204 0.0924 0.3588 
0.0291 0.0863 0.2221 
0.0656 0.1187 0.3289 


Total (3000+6000+9000+13500) 


0.3880 1.2115 2.7159 




0.1794 0.3552 0.5766 


Probability 


25% 50% 25% 



Table 8 Final turbulence integrals for the FA thin-layer model-OTP, J in units of 10 ^m 1 / 3 for SSO 
(Run 1-8: May 2005 to June 2006). 



Table 9 Bufton wind model. 





v G 


Vt 


Ht 


Lt 


Good 


2 


30 


9000 


4000 


Typical 


5 


35 


10000 


5000 


Bad 


8 


40 


11000 


6000 
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model-OTP layer heights for the ground-layer and free- 
atmosphere are used in the Button equation to obtain 
the model layer wind speed values. The model relative 
layer wind directions are obtained by using the model- 
OTP layer heights in the wind direction, 4>(h), defined 
in Equation 1 10| The wind speed model and wind di- 
rection model values for the SSO (Run f-8: M ay 2 005 
to June 2006) model-OTP are shown in Figure 13 and 
are tabulated in Table \W\ 



6.3 Summary Model-OTP 

The summary model-OTP (Kolmogorov) is the com- 
bination of the 'good', 'typical' and 'bad' GL and FA 
model-OTPs resulting in a set of nine possible thin- 
layer turbulence profiles. Each turbulence profile of 
the model-OTP includes information about the turbu- 
lence layer strength, wind speed and wind direction. 

A detailed summary table of the summary model- 
OTPs describing the observations at Siding Spring Ob- 
servatory, Australia for the SSO (Run 1-8: May 2005 



to June 2006) runs are listed in Table 11 (turbulence 
integrals); Table T2| (turbulence layer wind speeds) and 
Table 13 (turbulence layer wind direction). Note the 
model-OTP seeing, e, for the typical GL+FA of 1.24" 
is similar to historical median DIMM seeing of 1.25" 



reported by |Wood et al.| j |1995[ | 
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Figure 13 Wind model-OTP for the SSO (Run 1-8: May 2005 to June 2006) runs: squares 'good'; 
diamonds 'typical' and circles 'bad' (a) Button wind profile (b) Wind direction (empirical). Layer heights 
for model-OTP for SSO (Run 1-8: May 2005 to June 2006) model are marked as symbols. 







Button wind speed model (m/s) 


Wind direction 




height (m) 


good 


typical 


bad 


(degrees, ref 0m) 


GL 


37.5 


2.1981 


5.6605 


9.4200 


1.1674 




250 


2.2506 


5.7810 


9.6142 


7.5347 


FA 


1000 


2.5495 


6.3707 


10.4871 


26.9588 




3000 


5.1620 


9.9300 


14.7605 


61.3717 




6000 


19.0935 


23.4552 


27.9741 


85.7917 




9000 


32.0000 


38.6276 


43.7936 


95.5085 




13500 


10.4619 


26.4419 


41.6249 


100.3182 



Table 10 Tabulated values for the model wind speed and model wind direction for the GL and FA 
model-OTP layers for SSO (Run 1-8: May 2005 to June 2006) model. 



5: 

c 
G- 



Parameter 



37.5 

250 

1000 

3000 

6000 

9000 

13500 



Jgl 
Jfa 
J 



Fgl 
F fa 

(■FA 

e 



Probability 



GL 



FA 



Units 



l -13 m l/3 
10 -13 m l/3 
10 -13 m l/3 
10 -13 m l/3 
l -13 m l/3 
10 -13 m l/3 
10 -13 m l/3 



10 -13 m l/3 
l -13 m l/3 
10 -13 m l/3 



/1.0 

/i-o 

arcsecs 
arcsecs 
arcsecs 
arcsecs 
ms 



/1.0 



Model Turbulence Profiles (J, Integral) - SSO (Run 1-8: May 2005 to June 2006) 



Good 



Good 



1 



Typical 



Bad 



4.7129 
0.2485 
0.2413 
0.0315 
0.0204 
0.0291 
0.0656 



4.7129 
0.2485 
0.4772 
0.4368 
0.0924 
0.0863 
0.1187 



4.7129 
0.2485 
0.6677 
1.1384 
0.3588 
0.2221 
0.3289 



4.9614 
0.3880 
5.3493 



4.9614 
1.2115 
6.1729 



4.9614 
2.7159 
7.6773 



0.9275 
0.0725 
0.8277 
0.1794 
0.8659 
6.4233 
11.7922 



0.8037 
0.1963 
0.8277 
0.3552 
0.9436 
3.7172 
5.3516 



0.6462 
0.3538 
0.8277 
0.5766 
1.0755 
2.0123 
2.3291 



0.0625 0.1250 0.0625 



Typical 



Good 



Typical 



Bad 



6 



7.3989 
1.1574 
0.2413 
0.0315 
0.0204 
0.0291 
0.0656 



7.3989 
1.1574 
0.4772 
0.4368 
0.0924 
0.0863 
0.1187 



7.3989 
1.1574 
0.6677 
1.1384 
0.3588 
0.2221 
0.3289 



8.5562 
0.3880 
8.9442 



8.5562 
1.2115 
9.7678 



8.5562 
2.7159 
11.2721 



0.9566 
0.0434 
1.1478 
0.1794 
1.1787 
6.3684 
4.5855 



0.8760 
0.1240 
1.1478 
0.3552 
1.2427 
3.7043 
3.5038 



0.7591 
0.2409 
1.1478 
0.5766 
1.3542 
2.0098 
2.0310 



0.1250 0.2500 0.1250 



Bad 



Good 



7 



Typical 



9.4363 
3.8486 
0.2413 
0.0315 
0.0204 
0.0291 
0.0656 



9.4363 
3.8486 
0.4772 
0.4368 
0.0924 
0.0863 
0.1187 



13.2849 

0.3880 

13.6728 



13.2849 

1.2115 

14.4964 



0.9716 
0.0284 
1.4945 
0.1794 
1.5206 
6.2255 
2.2067 



0.9164 
0.0836 
1.4945 
0.3552 
1.5749 
3.6700 
2.0112 



Bad 



9 



9.4363 
3.8486 
0.6677 
1.1384 
0.3588 
0.2221 
0.3289 



13.2849 

2.7159 

16.0008 



0.8303 
0.1697 
1.4945 
0.5766 
1.6710 
2.0030 
1.5242 



0.0625 0.1250 0.0625 



Table 11 Tabulated values for the final model-OTP for SSO (Run 1-8: May 2005 to June 2006), with layers specified as turbulence integral, J, in units 



c 
e 

3 

CO 



lO-^mVs. 



to 
o 



Parameter 

37.5 

250 

1000 

3000 

6000 

9000 

13500 



Jgl 
Jfa 
J 



F GL 
F fa 



e 



70 



Probability 



GL 



FA 



Units 

m/s 

m/s 

m/s 

m/s 

m/s 

m/s 

m/s 



10 -13 TO l/3 

lO" 13 ™ 1 / 3 

l -13 m l/3 



/1.0 
/1.0 



arcsecs 
arcsecs 
arcsecs 



arcsecs 
ms 



/1.0 



Model Turbulence Profiles (Layer wind speeds) - 
Good Typical 



Good 



1 

2.1981 

2.2506 

2.5495 

5.1620 

19.0935 

32.0000 

10.4619 



Typical 



2 

2.1981 

2.2506 

6.3707 

9.9300 

23.4552 

38.6276 

26.4419 



Bad 



3 

2.1981 

2.2506 

10.4871 

14.7605 

27.9741 

43.7936 

41.6250 



4.9614 
0.3880 
5.3493 



4.9614 
1.2115 
6.1729 



4.9614 
2.7159 
7.6773 



0.9275 
0.0725 



0.8037 
0.1963 



0.6462 
0.3538 



0.8277 
0.1794 
0.8659 



0.8277 
0.3552 
0.9436 



0.8277 
0.5766 
1.0755 



6.4233 
11.7922 



3.7172 
5.3516 



2.0123 
2.3291 



0.0625 0.1250 0.0625 



SSO (Run 1-8: May 2005 to June 2006) 
Bad 



Good 



4 

5.6605 

5.7810 

2.5495 

5.1620 

19.0935 

32.0000 

10.4619 



Typical Bad 



Good 



5 

5.6605 

5.7810 

6.3707 

9.9300 

23.4552 

38.6276 

26.4419 



6 

5.6605 

5.7810 

10.4871 

14.7605 

27.9741 

43.7936 

41.6250 



7 

9.4200 

9.6142 

2.5495 

5.1620 

19.0935 

32.0000 

10.4619 



Typical 



9.4200 

9.6142 

6.3707 

9.9300 

23.4552 

38.6276 

26.4419 



8.5562 
0.3880 
8.9442 



8.5562 
1.2115 
9.7678 



8.5562 
2.7159 
11.2721 



13.2849 

0.3880 

13.6728 



13.2849 

1.2115 

14.4964 



0.9566 
0.0434 



0.8760 
0.1240 



0.7591 
0.2409 



0.9716 
0.0284 



0.9164 
0.0836 



1.1478 
0.1794 
1.1787 



1.1478 
0.3552 
1.2427 



1.1478 
0.5766 
1.3542 



1.4945 
0.1794 
1.5206 



1.4945 
0.3552 
1.5749 



6.3684 
4.5855 



3.7043 
3.5038 



2.0098 
2.0310 



6.2255 
2.2067 



3.6700 
2.0112 



0.1250 0.2500 0.1250 0.0625 0.1250 



Bad 



9 

9.4200 

9.6142 

10.4871 

14.7605 

27.9741 

43.7936 

41.6250 



13.2849 

2.7159 

16.0008 



0.8303 
0.1697 



1.4945 
0.5766 
1.6710 



2.0030 
1.5242 



0.0625 



Table 12 Tabulated values for the final model-OTP for the SSO (Run 1-8: May 2005 to June 2006), with layers specified as wind speeds. 



c 
2! 
o" 

53 
g 

8, 



o 
£. 

? 

CO" 

o 
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Model Turbulence Profiles (Layer wind directions) - SSO (Run 1-8: May 2005 to June 2006) 





GL 


Good 


Typical 


Bad 




FA 


Good 


Typical 


Bad 


Good 


Typical 


Bad 


Good 


Typical 


Bad 


Parameter 


Units 


1 


2 


3 


4 


5 


6 


7 


8 


9 


37.5 


degrees 


1.1674 


1.1674 


1.1674 


1.1674 


1.1674 


1.1674 


1.1674 


1.1674 


1.1674 


250 


degrees 


7.5347 


7.5347 


7.5347 


7.5347 


7.5347 


7.5347 


7.5347 


7.5347 


7.5347 


1000 


degrees 


26.9588 


26.9588 


26.9588 


26.9588 


26.9588 


26.9588 


26.9588 


26.9588 


26.9588 


3000 


degrees 


61.3717 


61.3717 


61.3717 


61.3717 


61.3717 


61.3717 


61.3717 


61.3717 


61.3717 


6000 


degrees 


85.7917 


85.7917 


85.7917 


85.7917 


85.7917 


85.7917 


85.7917 


85.7917 


85.7917 


9000 


degrees 


95.5085 


95.5085 


95.5085 


95.5085 


95.5085 


95.5085 


95.5085 


95.5085 


95.5085 


13500 


degrees 


100.3182 


100.3182 


100.3182 


100.3182 


100.3182 


100.3182 


100.3182 


100.3182 


100.3182 


Jgl 


10 -13 m l/3 


4.9614 


4.9614 


4.9614 


8.5562 


8.5562 


8.5562 


13.2849 


13.2849 


13.2849 


Jfa 


10~ 13 to 1/3 


0.3880 


1.2115 


2.7159 


0.3880 


1.2115 


2.7159 


0.3880 


1.2115 


2.7159 


J 


l -13 m l/3 


5.3493 


6.1729 


7.6773 


8.9442 


9.7678 


11.2721 


13.6728 


14.4964 


16.0008 


F GL 


/1.0 


0.9275 


0.8037 


0.6462 


0.9566 


0.8760 


0.7591 


0.9716 


0.9164 


0.8303 


Ffa 


/1.0 


0.0725 


0.1963 


0.3538 


0.0434 


0.1240 


0.2409 


0.0284 


0.0836 


0.1697 


zgl 


arcsecs 


0.8277 


0.8277 


0.8277 


1.1478 


1.1478 


1.1478 


1.4945 


1.4945 


1.4945 




arcsecs 


0.1794 


0.3552 


0.5766 


0.1794 


0.3552 


0.5766 


0.1794 


0.3552 


0.5766 


e 


arcsecs 


0.8659 


0.9436 


1.0755 


1.1787 


1.2427 


1.3542 


1.5206 


1.5749 


1.6710 


Oo 


arcsecs 


6.4233 


3.7172 


2.0123 


6.3684 


3.7043 


2.0098 


6.2255 


3.6700 


2.0030 


T 


ms 


11.7922 


5.3516 


2.3291 


4.5855 


3.5038 


2.0310 


2.2067 


2.0112 


1.5242 


Probability 


/1.0 


0.0625 


0.1250 


0.0625 


0.1250 


0.2500 


0.1250 


0.0625 


0.1250 


0.0625 



Table 13 Tabulated values for the final model-OTP for the SSO (Run 1-8: May 2005 to June 2006), with layers specified as wind directions. 
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7 Conclusions 



Acknowledgments 



This paper has reported on our turbulence profiling ob- 
servational results performed at SSO. A summary of 
results with example data has been presented for mea- 
surement spanning years 2005 to 2006 using the 7x7 
(Runs 1-6) and 17 x 17 (Runs 7-8) SLODAR instru- 
ments. The observational results has facilitated the 
site-characterisation of the optical turbulence profile 
with the implementation of model-OTP (Kolmogorov) 
for SSO. The model-OTP describes the turbulence layer 
strength (Table |TTj ), layer wind speed (Table [l2| and 
layer wind direction (Table 13 1. The model-OTP is 
useful for the prediction of adaptive optics performance 
at SSO (forthcoming paper) using simulation codes. 
Prior to the commencement of our SLODAR campaign 
the seeing statistics were relatively well understood 
with DIMM seeing measurements ( |Wood et al.| 1995 1 , 
whereas the vertical structure of the turbulence pro- 
file structure was relatively unknown, with only a few 
SCIDAR profiles reported by|Klueckers et al.| (fl998 >. 



The following conclusions can be stated about the at- 
mospheric turbulence above SSO based on our data: 

• A measured median atmospheric seeing of around 
1.2" (Kolmogorov model with mirror/dome see- 
ing removed). The seeing is discussed in Sec- 
tion o 

• Ground-layer turbulence dominates, with ~80% 
turbulence below 500 m (h 5 oo). The structure 
of the turbulence is discussed in Section 15.11 
The presence of a strong ground-layer was also 
discovered during site measu rements at Moun t 
John University Observatory ( Mohr et al.|2010 1 . 
The dominate ground-layer has promising im- 
plications for GLAO. A forthcoming paper will 
investigate the performance of GLAO based on 
the model-OTP for SSO. 

• Non-Kolmogorov turbulence is observed espe- 
cially for ground-layer, with /3 avg ~ 10/3. The 
power-law slope of the spatial phase fluctuations 
is discussed in Section |5TT] The non-Kolomogorov 
spectrum results in a different scaling of seeing 
with wavelength to that conventionally assumed. 

• Turbulence profiles shows a number of dynam- 
ical characteristics; (i) most intense near the 
ground (below 100 m); (ii) fluctuate in intensity, 
appearing in 'bursts' with timescales of several 
minutes. These characteristics are shown in Fig- 
ure [6] as a sequence of turbulence profiles. 

• Mirror/dome seeing can be a significant fraction 
of the ground-layer turbulence. This is evident 
as a zero-height static contribution in the tur- 
blent l ayer wind spee d m easurements (refer to 
Section 



5.2 and Figure 10 1 



The free-atmosphere turbulence is comparable 
to 'good' seeing sites. From the model-OTP 
shown in Table |11| the free-atmosphere seeing 
for SSO is 0.18 (25%), 0.36 (50%) and 0.58 (75%) 
arcsecs. The Cerro Pachon (seeing^0.75 ) mode l 
(CP Model, |Tokovinin and Travouillonl ( |2006[ )) 
for the free-atmosphere reports 0.29 (25%), 0.40 
(50%) and 0.55 (75%) arcsecs. 
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mitted 2009) conducted at the Research School of As- 
tronomy and Astrophysics (RSAA) of the Australian 
National University (ANU). We thank Peter Conroy 
(ANU) for the design and manufacture of the instru- 
mentation. The authors are grateful to Jon Lawrence 
of the Australian Astronomical Observatory ( AAO) for 
his suggestions regarding the original version of the 
manuscript and to the referees for their comments. 
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